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ON ESTIMATING THE VENUS SPIN VECTOR FROM 
DATA OBTAINED DURING THE PLANETARY 
EXPLORER MISSION 


Peter Argentlero 
George Wyatt 

SUMMARY 

This paper demonstrates the feasibility of obtaining improved estimates of 
the Venus Spin Vector through a least squares processing of range rate data 
from P.E. probes on the surface of Venus. It is shown that the probes need 
transmit data only for a few moments after landing in order to obtain accurate 
estimates. This procedure estimates the right ascension of the spin vector 
better than the declination. The estimation procedure remains viable even if a 
few of the probes do not survive impact. Some improvement in the estimates of 
probe locations is also obtainable. 
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ON ESTIMATING THE VENUS SPIN VECTOR FROM 


DATA OBTAINED DURING THE PLANETARY 
EXPLORER MISSION 


INTRODUCTION 

The Planetary Explorer (P.E.) mission is designed to land a total of eight 
probes on the surface of Venus during the 1977 Venus opportunity. All these 
probes will be equipped with devices for returning range rate data to the Earth. 
These devices are expected to survive impact and to continue transmitting data 
from the planet's surface. The Venus spin vector as well as the positions of the 
probes on the surface are observable in this data. Thus the correct processing 
of this data should provide at least some improvement in our knowledge of these 
parameters. This prospect is an attractive one because of the considerable 
interest in the spin vector of Venus and because the data necessary for the esti- 
mation is obtained as a byproduct of the P.E. mission and thus its acquisition 
imposes no further constraints or compromises on mission planning. 

The purpose of this paper is to indicate the sort of improvement in spin vec- 
tor and probe position estimates one may reasonably expect from the processing 
of such data. This was done by duplicating the ensemble calculations associated 
with a weighted least squares with a-priori estimation technique applied to range 
rate data which was assumed to be unbiased and uncorrelated. The weighting 
matrix was assumed to be the inverse of the covariance matrix of the noise on 
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the data. Attention is focused primarily on the spin vector estimation. After a 
brief discussion of previous efforts to estimate the Venus spin vector, a presen- 
tation and analysis of the Authors* results will be given. 

PREVIOUS ESTIMATES OF VENUS SPIN VECTOR 

Several radar determinations of the Venus spin vector have been reported 
in the literature. This technique involves the transmission of a cw signal to the 
planet and the analysis of the power spectrum of the return signal. If the trans- 
mitted frequency corrected for relative motion between Venus and the tracking 
station is used as the zero point of the spectrum, then every other frequency can 
be interpreted as a doppler shift, and the intensity associated with that frequency 
can be related to a line of constant radial velocity on the planet. These lines are 
easy to describe. A plane going through the tracking station and parallel to the 
spin axis of Venus and intersecting the planet*s surface describes a circle. That 
segment of the circle visible to the tracking station is a member of the family of 
lines of constant radial velocity. A frequency at which a peak in the spectra oc- 
curs corresponds to a region on Venus which is rougher than adjacent regions. 
Thus a cw spectrum may be viewed as a type of map of the radar brightness of 
the surface of the planet. By obtaining such spectra at different times it should 
be possible to track specific surface features and from such information to infer 
the spin vector of the planet under investigation. Carpenter [1] reproduces sev- 
eral of these radar spectra and identifies several surface features as peaks in 
the spectra. The mathematics of how one might utilize such information in order 
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to estimate a spin vector is given in some detail by Shapiro [2]. Basically a 
weighted least squares estimation procedure is used with the surface features 
associated with peaks in the spectra treated as point sources of range rate data. 
This method has been applied by Carpenter [3], and [4] , Smith [5]. Goldstein [6] 
and Shapiro [2]. Their various estimates of the rotation period of Venus are 
quite close to the so-called synodic resonance period of 243.16 days. If Venus 
were to have this period, it would rotate backward relative to the Earth four 
times between each inferior conjunction and thus present the same side of the 
Earth at each inferior conjunction. The uncertainties associated with present 
estimates are such that the hypothesis that Venus is in synodic resonance is 
still unsettled. 

Another interesting aspect of the radar determinations of the Venus spin 
vector is that statistically speaking they are mutually incompatible. Table 1 
exhibits spin vector estimates and associated statistics obtained from four 


Table 1 

Recent Radar Determinations and Associated Standard 
Deviations of Venus Spin Vector 


Determination 

Rotation Period 
(In Days) 

Right Ascension 
in Degrees 

Declination 
in Degrees 

1. Carpenter (1970 

242.98 ± .04 

94 ±3 

-71.5 ± 1 

2. Shapiro (1967) 

243.09 ± .18 

84.7 ± 1.8 

-65.8 ± 1.2 

3. Dyce, Pettengill, 




and Shapiro (1967) 

244.3 ± 2 

90.9 ± 1 

-66.4 ± 1 

4. Goldstein (1967) 

242.6 ± .6 

98 ±5 

-69 ± 2 
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different radar determinations; Carpenter [l], Shapiro [7], Dyce, Pettengill, and 
Shapiro [8], and Goldstein [9]. 

Clearly the standard deviations of Table 1 give a too optimistic representa- 
tion of the accuracy of the radar determination technique. The reason is that 
ensemble calculations giving standard deviations of the least squares estimation 
process do not reflect the impact of modeling errors. The most obvious modeling 
error associated with the technique is that of treating the source of a peak in a 
frequency spectrum as a point. Any surface feature on Venus capable of causing 
a significant perturbation in the return spectrum of a radar scan must be quite 
extensive. A more subtle though possibly less important modeling error is the 
assumption that the return rate data is uncorrelated. Corrections must be made 
to the doppler data before it can be used for least squares estimation. These 
corrections are responsible for a portion of the noise on the data. If this portion 
is significant, then it is improper to model the noise on the data as an uncorre- 
lated or "white” random process. The correct procedure would then be to solve 
for a bias and perhaps a scale factor error in the data. If this were done the 
ensemble calculation would yield a more accurate though quantitatively less im- 
pressive figure for the standard deviation of the Venus spin vector estimate. An- 
other possibility is that the cause of the peaks on the return spectra may be a 
complex combination of phenomena rather than a fixed surface feature. 

It is not easy to obtain reasonable estimates of the quality of the radar 
determination technique. Carpenter [1] suggests that the formal standard deviation 
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numbers associated with radar determinations should be increased several times. 


This is quite vague, but it would be imprudent to be more precise. 

SPIN VECTOR ESTIMATION WITH P.E. PROBES 

The P.E. mission will consist of two separate launchings of multiprobe 
spacecraft. Each launch will land a main probe equipped with a two way doppler 
devices and three mini probes each equipped with a one way doppler device. The 
geometric distribution of the probes is given by Figure 2. The succeeding 
analysis rests on these modeling assumptions. 

(1) For the duration of their transmission of doppler data, the probes re- 
main stationary relative to the Venus surface. 

(2) For this same duration the Venus spin vector is constant. 

(3) The noise on the range rate data transmitted by the probes is stationary 
and uncorrelated. 

(4) Unbiased estimates of the Venus spin vector and the locations and ef- 
fective radii of the probes are available. It is also assumed that these 
estimates are statistically uncorrelated and that their standard devia- 
tions are known. 

Assumptions one and two should disturb no one. Assumption three is some- 
what more troublesome. We defer discussion of this assumption to a later stage 
of the analysis. Concerning assumption four, the estimate of the Venus spin 
vector and its associated uncertainty would presumably be borrowed from a 
radar determination. The probe positions and associated uncertainties would be 
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obtained from extensive post-flight analysis of all relevant data gathered during 
the P.E. mission. This would include accelerometer, temperature and pressure 
measurements and range rate data. 

The conventional way to process the doppler data from the P.E. probes is 
to form the usual weighted least squares with a-priori loss function and choose 
the spin vector and probe positions which minimize this loss function. To be 
specific, let X be a vector of dimension twenty seven and whose elements are 
estimates of the spin vector and positions of each of the eight probes. Let the 
vector Y be the range rate measurements, arranged in some sequence, which 
would be obtained in the absence of noise if X contained the true values of the 
parameters in question. The vector Y is a known function of X symbolized oy 

Y = f (X) (1) 

Let Q be the covariance matrix of the noise on the observations and let P be the 
covariance matrix of an a-prori estimate X’ of the parameters. Then the loss 
function is 


L(X) = (V - f (X)) T Q ■* (Y - f (X)) + (X' - X) 5 P' 1 (X' - X) (2) 

where Y is the vector of measured values of range rate data obtained from the 
probes. Notice that by assumptions three and four the matrices Q and P are 
diagonal - a fact which greatly simplifies following computations. The least 
squares estimate is defined as that value of X which minimizes the loss function 
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of Equate 2. Since the data necessary to implement this procedure is not yet 
available, interest is focused on just the statistical properties of this estimation 
procedure. The covariance matrix of the least squares estimator may be ob- 
tained if one more assumption is imposed. It must be supposed that the function 
of equation 1 can be accurately represented as a first ordei Taylor series about 

A A 

X where X is the least squares estimate. Thus we assume that equation 1 can 
be written as 


Y - f (X) = A (X - X) (3) 

The symbol A represents the variational matrix. If N be the total number of 
range rate measurements, then A is an N by 27 dimensional matrix. The ele- 
ment in the i-th row and j-th column of A is the partial derivative of the i-th 

\ 

component of Y with respect to the j-th component of X. It ie relatively easy to 
obtain an analytic?! expression for A. The details are found in appendix one. If 
equation 3 is valid, then the covariance matrix of the least squares estimator is 
given by 

COV (X) = (A T Q- 1 A + p- 1 ) -1 (4) 

Equation 4 provides a mode for the performance of a parametric study of 
the accuracy attainable in the estimation of the Venus spin vector. Interest was 
focused primarily on the variation of this accuracy with respect to variations in 
the following parameters 
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(1) Length of time the probes survive on the surface 

(2) Quality of a-prori information 

(3) Size of the noise on the data 

This paper reports on the results of such a parametric study* As usual, nomi- 
nal values were established for all relevant parameters. Certain parameters 
were then systematically varied with the other parameters fixed at their 
nominal values. Equation 4 is then used to obtain the corresponding standard 
deviations of the spin vector and probe locations. There is a certain amount of 
arbitrariness involved in the selection of nominal values. It is not easy, for 
instance, to decide w'hat are reasonable values for the a-priori uncertainties of 
the spin vector at somo time several years in the future when the P.E. mission 
is to be executed. Also the time at which the probes land is a factor since this 
affects tracking sight geometry. This too is somewhat arbitrary although one 
would of course choose a time during which dual coverage from the D.S.N. net- 
work is possible. The data chosen for the landing is May 19, 1977. Dual cover- 
age from Goldstone and Madrid is assumed. The nominal value of the spin vec- 
tor as obtained from Cl] is 242.982 days for the period and 94.1 and -71.4 degrees 
respectively for the right ascension and declination of the spin vector (equator of 
1950). The longitude and latitude of each probe is given in target coordinates in 
[10]. The data acquisition rate is assumed to be one per minute. It is also as- 
sume that all biases are estimated in the least squares procedure. 

A summary of a-priori standard deviations is provided in Table II. 
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Table 2 

A Priori Standard Deviations 


Parameter 

Standard Deviation 

Right ascension of main probes 

.115° 

Declination of main probes 

.115° 

Right ascension of mini probes 

.165° 

Declination of mini probes 

.165° 

Range rate from main probes 

5 mm/sec. 

Range rate from mini probes 

5 cm/sec. 

Period of Venus rotation 

1 day 

Right ascension of Venus spin vector 

10° 

Declination of Venus spin vector 

5° 

Effective radius of Venus at each probe 

5 km 


The usual standard deviation figure for two way range rate and for a one 
per minute sampling rate is 1 mm/sec. See for instance Blackshear and 
Williams [11] who use this figure for two way range rate data in an error study 
similar to this one but for the Viking project. The larger figure of 5 mm/sec. 
for the two way range rate data on the main probes was used to compensate for 
the fact that the data was assumed to have no time correlation. Since the vari- 
ous corrections which must be made to the data in fact do tend to introduce time 
correlations this assumption is a questionable one. Its legitimacy can only be 
defended if a bias on the data is also solved for in the estimation procedure. In 
this error study no provisions have been made for bias estimation. Hence it was 
felt that to use the usual standard deviation figure for the noise on the two way 
range rate data would lead to an excessively optimistic result. The use of a 
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larger standard deviation number was an effort to compensate for this optimism. 
Since the presence of biases on one way doppler data is a much more serious 
problem, the standard deviation of the noise on the mini probe data was set at 
ten times the corresponding figure for the noise on the main probe data. 

Figures 2, 3 and 4 display the possibilities of estimating the period, right 
ascension, and declination of the Venus spin vector as a function of the length 
of time the probe's range rate devices survive on the surface. Since it is not 
certain that all range rate devices will survive parachute openings or impact on 
the surface, the figures also demonstrate the deterioration in the estimation 
procedure if a mini probe per launch fails to function and also if a mini probe 
and a main probe on each launch fail to function. A salient feature of the figures 
is that most of the improvement in the estimates occurs within the first few 
minutes after impact. This fact suggests that the primary reason for the feasi- 
bility of this estimation procedure is the coupling provided by Equation 4 between 
the spin vector estimate and the relatively small standard deviations of the 
a priori estimates of probe locations. If this were true one would expect the 
quality of the estimation procedure to be far more sensitive to changes in the 
standard deviations of the a priori estimates of probe locations than to changes 
in the standard deviations of the noise on the data. This appears to be the case. 
For example, with regard to the nominal case with all probes lasting thirty 
minutes the standard deviation of the rotation period is .28 days. When the 
standard deviations of the noise on the data are doubled the number becomes .31 
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days. With the standard deviations of the noise returned to nominal values but 


with tL> standard deviations of the a priori probe location estimates doubled the 
standard deviation of the estimate of the period rises to .4 days. 

It is apparent from a comparison of figures 3 and 4 that the right ascension 
of the Venus spin vector is far more observable in this experiment than the 
declination. A glance at Table 1 rev sals that the situation is precisely the op- 
posite with regard to the radar determination procedure. It is not obvious why 
this should be so. But in this sense, at least, the two estimation procedures 
should neatly compliment each other. 

The ability of this estimation procedure to improve knowledge of probe posi- 
tions is somewhat less impressive. After thirty minutes of tracking, knowledge 
of the positions of the main probes is improved by appi imately 15%. Knowl- 
edge of the positions of the mini probes after the same period of tracking is 
improved by approximately 10%. After an hour of tracking these percentage 
improvements are respectively 18% and 15%. 

CONCLUSION 

This paper has demonstrated the feasibility of utilizing range rate data 
generated by P.E. probes in order to estimate the Venus spin vector. The esti- 
mation procedure is viable even if each probe transmits doppler data for just a 
few minutes after impact. The standard deviation figures associated with this 
estimation procedure are not dependent on questionable modeling assumptions 
for their validity. Consequently they are a true measure of the estimated accuracy. 
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The right ascension of the spin vector is far more observable in the data 
than the declination. Fortuitously, the radar determination technique estimates 
the declination better than the right ascension. Thus the two estimation tech- 
niques compliment each other. 

The improvements in probe position uncertainty are somewhat marginal. 

Finally, one of the happy aspects of this proposed experiment is that it 
utilizes data which is generated essentially as a byproduct of the P.E. mission 
and hence its implementation poses no additional constraints on P.E. mission 
planning. 
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Figure 1. Position of Probes on Venus Surface 
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Figure 2. Standard Deviation of Period as Function of Time 
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Figure 3. Standard Deviation of Right Ascension as a Function of Time 
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Figure 4. Standard Deviation of Declination os o Function of Time 
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APPENDIX 

THE VARIATIONAL MATRIX FOR SPIN VECTOR ESTIMATION 

Let Y be a vector of range rate measurements obtained from P.E. probes 
on the surface of Venus. Let X be the set of parameters that are to be estimated. 
In this case X consists of the three components of the spin vector and three 
components for each of the probes on the surface. The variational matrix A is 
defined by specifying that the element in the i-th row and j-th column of A is the 
partial derivative, evaluated at nominal values of X, of the i-th component of Y 
with regard to the j-th component of X. Notice that a given range rate measure- 
ment has a non- zero partial derivative only with regard to the components of the 
spin vector and the position components of the probe which generated the meas- 
urement. What follows is a derivation of these six partial derivatives of a range 
rate measurement from any given probe and at any given time. 

Let R q be the position vector of a probe at time of impact t Q . If R(t) repre- 
sents the position of the probe at a later time t and if w represents the spin vec- 
tor of Venus, then the setor equation of motion of the probe can be written as 

R(t) =n^R(t) (1) 

where 
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The unique solution to this matrix differential equation is 


R(t) = e n ' < ‘' , » ) R. 


( 2 ) 


where by definition 


Q w (t-t 0 ) = j 


N 

Lim y 


k! 


(3) 


The right side of Equation 3 can be written in the following closed form 




e ^i ( t-t o ) 


(n w -.x 2 i) (fl w -\ a i> 

(\j - \ a ) - x 3 > 


(4) 


+ e 




(\ 2 -\ 3 ) +e (x 3 -^,) (S-\ 2 ) 


where \ 1# k 2 , X 3 are the eigen values of fl w . These eigen values are known to 
be 0, i|w|, - i|w|. Hence after some manipulation, Equation 4 can be written as 


e 





“■ cos |w| (t - t 0 ) + sin |w| (t - t 0 ) 

|w| 


By differentiating Equation 2 one obtains 


(5) 


R(t) =fi„ e 51 * 0 *'*’ R 0 . 


( 6 ) 


Letting the symbol B t represent the unit vector from the position of the probe 
at time t to the tracking station, the scaler range rate D(t) measured at the 
tracking station at time t is 
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( 7 ) 


D(t) = B* R(t) = B* fi. e n * (, ‘*o ) R 0 . 

Thus the gradient of D(t) with regard to R 0 is 

V >0 D(t)=B{ fi, (8) 

To obtain the gradient of D(t) with respect to W, one begins by writing Equation 
1 in the following form 

R (‘> = n S(t) w < 9 > 

where 

/o -R,(t) R y (t)\ 

n R ( „= R .(‘) 0 

\-S y (t) R„(t) 0 J 


Hence 

D(t) = Bjn„ (t) w (1.0) 

As long as t “ t 0 does not represent a duration of more than a few hours, the 
dependence of the components of n R on W can be neglected and the gradient of 
D(t) with regard to W can be stated as 

V w D(t)=B*n B<t) . (ID 

To see that this is true consider the correct equation for the partial derivative 


of D(t) with respect to W x 
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It is easy to show that the first two terms on the right side of Equation 12 
are the only ones that need be considered. Consider for instance the third term 
on the right side of Equation 12. The nominal value of W x is approximately 
3 (10)’ 4 radians/hr. The terms 


3R 


3R 

z 

and 

y 

3w , 

3w 

X 


are bounded respectively by |r J (t - t 0 ) and I R y I (t - t 0 ) where t - to is meas- 
ured in hours. Thus the first two terms are seen to be four orders of magnitude 
greater than succeeding terms. By deleting these smaller terms and performing 
similar deletions in the equations for the partial derivatives of D(t) with respect 
to W y and W z , one obtains the representation of Equation 11. It is interesting to 
observe that if the period of the Venus rotation were substantially shorter or if 
data from the P.E. probes were to be processed over a longer period of time, 
say several days, then that portion of the variational matrix obtained by Equation 
11 would be substantially more complicated. 
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